Quantification of mediation effects of white matter functional characteristics on cognitive decline in aging.

Cognitive decline with aging involves multifactorial processes, including changes in brain structure and function. This study focuses on the role of white matter functional characteristics, as reflected in blood oxygenation level-dependent signals, in age-related cognitive deterioration. Building on previous research confirming the reproducibility and age-dependence of blood oxygenation level-dependent signals acquired via functional magnetic resonance imaging, we here employ mediation analysis to test if aging affects cognition through white matter blood oxygenation level-dependent signal changes, impacting various cognitive domains and specific white matter regions. We used independent component analysis of resting-state blood oxygenation level-dependent signals to segment white matter into coherent hubs, offering a data-driven view of white matter's functional architecture. Through correlation analysis, we constructed a graph network and derived metrics to quantitatively assess regional functional properties based on resting-state blood oxygenation level-dependent fluctuations. Our analysis identified significant mediators in the age-cognition relationship, indicating that aging differentially influences cognitive functions by altering the functional characteristics of distinct white matter regions. These findings enhance our understanding of the neurobiological basis of cognitive aging, highlighting the critical role of white matter in maintaining cognitive integrity and proposing new approaches to assess interventions targeting cognitive decline in older populations.

Optimization of the flip angles of narrow-band editing pulses in J-difference edited MRS of lactate at 3T.

PURPOSE: Application of highly selective editing RF pulses provides a means of minimizing co-editing of contaminants in J-difference MRS (MEGA), but it causes reduction in editing yield. We examined the flip angles (FAs) of narrow-band editing pulses to maximize the lactate edited signal with minimal co-editing of threonine. METHODS: The effect of editing-pulse FA on the editing performance was examined, with numerical and phantom analyses, for bandwidths of 17.6-300 Hz in MEGA-PRESS editing of lactate at 3T. The FA and envelope of 46 ms Gaussian editing pulses were tailored to maximize the lactate edited signal at 1.3 ppm and minimize co-editing of threonine. The optimized editing-pulse FA MEGA scheme was tested in brain tumor patients. RESULTS: Simulation and phantom data indicated that the optimum FA of MEGA editing pulses is progressively larger than 180° as the editing-pulse bandwidth decreases. For 46 ms long 17.6 Hz bandwidth Gaussian pulses and other given sequence parameters, the lactate edited signal was maximum at the first and second editing-pulse FAs of 241° and 249°, respectively. The edit-on and difference-edited lactate peak areas of the optimized FA MEGA were greater by 43% and 25% compared to the 180°-FA MEGA, respectively. In-vivo data confirmed the simulation and phantom results. The lesions of the brain tumor patients showed elevated lactate and physiological levels of threonine. CONCLUSION: The lactate MEGA editing yield is significantly increased with editing-pulse FA much larger than 180° when the editing-pulse bandwidth is comparable to the lactate quartet frequency width.

Whole-brain, gray, and white matter time-locked functional signal changes with simple tasks and model-free analysis.

Recent studies have revealed the production of time-locked blood oxygenation level-dependent (BOLD) functional MRI (fMRI) signals throughout the entire brain in response to tasks, challenging the existence of sparse and localized brain functions and highlighting the pervasiveness of potential false negative fMRI findings. "Whole-brain" actually refers to gray matter, the only tissue traditionally studied with fMRI. However, several reports have demonstrated reliable detection of BOLD signals in white matter, which have previously been largely ignored. Using simple tasks and analyses, we demonstrate BOLD signal changes across the whole brain, in both white and gray matters, in similar manner to previous reports of whole brain studies. We investigated whether white matter displays time-locked BOLD signals across multiple structural pathways in response to a stimulus in a similar manner to the cortex. We find that both white and gray matter show time-locked activations across the whole brain, with a majority of both tissue types showing statistically significant signal changes for all task stimuli investigated. We observed a wide range of signal responses to tasks, with different regions showing different BOLD signal changes to the same task. Moreover, we find that each region may display different BOLD responses to different stimuli. Overall, we present compelling evidence that, just like all gray matter, essentially all white matter in the brain shows time-locked BOLD signal changes in response to multiple stimuli, challenging the idea of sparse functional localization and the prevailing wisdom of treating white matter BOLD signals as artifacts to be removed.

Spatiospectral image processing workflow considerations for advanced MR spectroscopy of the brain.

Magnetic resonance spectroscopy (MRS) is one of the few non-invasive imaging modalities capable of making neurochemical and metabolic measurements in vivo. Traditionally, the clinical utility of MRS has been narrow. The most common use has been the "single-voxel spectroscopy" variant to discern the presence of a lactate peak in the spectra in one location in the brain, typically to evaluate for ischemia in neonates. Thus, the reduction of rich spectral data to a binary variable has not classically necessitated much signal processing. However, scanners have become more powerful and MRS sequences more advanced, increasing data complexity and adding 2 to 3 spatial dimensions in addition to the spectral one. The result is a spatially- and spectrally-variant MRS image ripe for image processing innovation. Despite this potential, the logistics for robustly accessing and manipulating MRS data across different scanners, data formats, and software standards remain unclear. Thus, as research into MRS advances, there is a clear need to better characterize its image processing considerations to facilitate innovation from scientists and engineers. Building on established neuroimaging standards, we describe a framework for manipulating these images that generalizes to the voxel, spectral, and metabolite level across space and multiple imaging sites while integrating with LCModel, a widely used quantitative MRS peak-fitting platform. In doing so, we provide examples to demonstrate the advantages of such a workflow in relation to recent publications and with new data. Overall, we hope our characterizations will lower the barrier of entry to MRS processing for neuroimaging researchers.

Automatic Preprocessing Pipeline for White Matter Functional Analyses of Large-Scale Databases.

Recently, increasing evidence suggests that fMRI signals in white matter (WM), conventionally ignored as nuisance, are robustly detectable using appropriate processing methods and are related to neural activity, while changes in WM with aging and degeneration are also well documented. These findings suggest variations in patterns of BOLD signals in WM should be investigated. However, existing fMRI analysis tools, which were designed for processing gray matter signals, are not well suited for large-scale processing of WM signals in fMRI data. We developed an automatic pipeline for high-performance preprocessing of fMRI images with emphasis on quantifying changes in BOLD signals in WM in an aging population. At the image processing level, the pipeline integrated existing software modules with fine parameter tunings and modifications to better extract weaker WM signals. The preprocessing results primarily included whole-brain time-courses, functional connectivity, maps and tissue masks in a common space. At the job execution level, this pipeline exploited a local XNAT to store datasets and results, while using DAX tool to automatic distribute batch jobs that run on high-performance computing clusters. Through the pipeline, 5,034 fMRI/T1 scans were preprocessed. The intraclass correlation coefficient (ICC) of test-retest experiment based on the preprocessed data is 0.52 - 0.86 (N=1000), indicating a high reliability of our pipeline, comparable to previously reported ICC in gray matter experiments. This preprocessing pipeline highly facilitates our future analyses on WM functional alterations in aging and may be of benefit to a larger community interested in WM fMRI studies.

Functional alterations in bipartite network of white and grey matters during aging.

The effects of normal aging on functional connectivity (FC) within various brain networks of gray matter (GM) have been well-documented. However, the age effects on the networks of FC between white matter (WM) and GM, namely WM-GM FC, remains unclear. Evaluating crucial properties, such as global efficiency (GE), for a WM-GM FC network poses a challenge due to the absence of closed triangle paths which are essential for assessing network properties in traditional graph models. In this study, we propose a bipartite graph model to characterize the WM-GM FC network and quantify these challenging network properties. Leveraging this model, we assessed the WM-GM FC network properties at multiple scales across 1,462 cognitively normal subjects aged 22-96 years from three repositories (ADNI, BLSA and OASIS-3) and investigated the age effects on these properties throughout adulthood and during late adulthood (age ≥70 years). Our findings reveal that (1) heterogeneous alterations occurred in region-specific WM-GM FC over the adulthood and decline predominated during late adulthood; (2) the FC density of WM bundles engaged in memory, executive function and processing speed declined with age over adulthood, particularly in later years; and (3) the GE of attention, default, somatomotor, frontoparietal and limbic networks reduced with age over adulthood, and GE of visual network declined during late adulthood. These findings provide unpresented insights into multi-scale alterations in networks of WM-GM functional synchronizations during normal aging. Furthermore, our bipartite graph model offers an extendable framework for quantifying WM-engaged networks, which may contribute to a wide range of neuroscience research.

Changes in white matter functional networks across late adulthood.

Introduction: The aging brain is characterized by decreases in not only neuronal density but also reductions in myelinated white matter (WM) fibers that provide the essential foundation for communication between cortical regions. Age-related degeneration of WM has been previously characterized by histopathology as well as T2 FLAIR and diffusion MRI. Recent studies have consistently shown that BOLD (blood oxygenation level dependent) effects in WM are robustly detectable, are modulated by neural activities, and thus represent a complementary window into the functional organization of the brain. However, there have been no previous systematic studies of whether or how WM BOLD signals vary with normal aging. We therefore performed a comprehensive quantification of WM BOLD signals across scales to evaluate their potential as indicators of functional changes that arise with aging. Methods: By using spatial independent component analysis (ICA) of BOLD signals acquired in a resting state, WM voxels were grouped into spatially distinct functional units. The functional connectivities (FCs) within and among those units were measured and their relationships with aging were assessed. On a larger spatial scale, a graph was reconstructed based on the pair-wise connectivities among units, modeling the WM as a complex network and producing a set of graph-theoretical metrics. Results: The spectral powers that reflect the intensities of BOLD signals were found to be significantly affected by aging across more than half of the WM units. The functional connectivities (FCs) within and among those units were found to decrease significantly with aging. We observed a widespread reduction of graph-theoretical metrics, suggesting a decrease in the ability to exchange information between remote WM regions with aging. Discussion: Our findings converge to support the notion that WM BOLD signals in specific regions, and their interactions with other regions, have the potential to serve as imaging markers of aging.

Experimental demonstration of diffusion limitations on resolution and SNR in MR microscopy.

PURPOSE: MR microscopy is in principle capable of producing images at cellular resolution (<10 µm), but various factors limit the quality achieved in practice. A recognized limit on the signal to noise ratio and spatial resolution is the dephasing of transverse magnetization caused by diffusion of spins in strong gradients. Such effects may be reduced by using phase encoding instead of frequency encoding read-out gradients. However, experimental demonstration of the quantitative benefits of phase encoding are lacking, and the exact conditions in which it is preferred are not clearly established. We quantify the conditions where phase encoding outperforms a readout gradient with emphasis on the detrimental effects of diffusion on SNR and resolution. METHODS: A 15.2 T Bruker MRI scanner, with 1 T/m gradients, and micro solenoid RF coils < 1 mm in diameter, were used to quantify diffusion effects on resolution and the signal to noise ratio of frequency and phase encoded acquisitions. Frequency and phase encoding's spatial resolution and SNR per square root time were calculated and measured for images at the diffusion limited resolution. The point spread function was calculated and measured for phase and frequency encoding using additional constant time phase gradients with voxels 3-15 µm in dimension. RESULTS: The effect of diffusion during the readout gradient on SNR was experimentally demonstrated. The achieved resolutions of frequency and phase encoded acquisitions were measured via the point-spread-function and shown to be lower than the nominal resolution. SNR per square root time and actual resolution were calculated for a wide range of maximum gradient amplitudes, diffusion coefficients, and relaxation properties. The results provide a practical guide on how to choose between phase encoding and a conventional readout. Images of excised rat spinal cord at 10 µm × 10 µm in-plane resolution demonstrate phase encoding's benefits in the form of higher measured resolution and higher SNR than the same image acquired with a conventional readout. CONCLUSION: We provide guidelines to determine the extent to which phase encoding outperforms frequency encoding in SNR and resolution given a wide range of voxel sizes, sample, and hardware properties.

Superficial white matter across development, young adulthood, and aging: volume, thickness, and relationship with cortical features.

Superficial white matter (SWM) represents a significantly understudied part of the human brain, despite comprising a large portion of brain volume and making up a majority of cortico-cortical white matter connections. Using multiple, high-quality datasets with large sample sizes (N = 2421, age range 5-100) in combination with methodological advances in tractography, we quantified features of SWM volume and thickness across the brain and across development, young adulthood, and aging. We had four primary aims: (1) characterize SWM thickness across brain regions (2) describe associations between SWM volume and age (3) describe associations between SWM thickness and age, and (4) quantify relationships between SWM thickness and cortical features. Our main findings are that (1) SWM thickness varies across the brain, with patterns robust across individuals and across the population at the region-level and vertex-level; (2) SWM volume shows unique volumetric trajectories with age that are distinct from gray matter and other white matter trajectories; (3) SWM thickness shows nonlinear cross-sectional changes across the lifespan that vary across regions; and (4) SWM thickness is associated with features of cortical thickness and curvature. For the first time, we show that SWM volume follows a similar trend as overall white matter volume, peaking at a similar time in adolescence, leveling off throughout adulthood, and decreasing with age thereafter. Notably, the relative fraction of total brain volume of SWM continuously increases with age, and consequently takes up a larger proportion of total white matter volume, unlike the other tissue types that decrease with respect to total brain volume. This study represents the first characterization of SWM features across the large portion of the lifespan and provides the background for characterizing normal aging and insight into the mechanisms associated with SWM development and decline.

Short superficial white matter and aging: a longitudinal multi-site study of 1293 subjects and 2711 sessions.

It is estimated that short association fibers running immediately beneath the cortex may make up as much as 60% of the total white matter volume. However, these have been understudied relative to the long-range association, projection, and commissural fibers of the brain. This is largely because of limitations of diffusion MRI fiber tractography, which is the primary methodology used to non-invasively study the white matter connections. Inspired by recent anatomical considerations and methodological improvements in superficial white matter (SWM) tractography, we aim to characterize changes in these fiber systems in cognitively normal aging, which provide insight into the biological foundation of age-related cognitive changes, and a better understanding of how age-related pathology differs from healthy aging. To do this, we used three large, longitudinal and cross-sectional datasets (N = 1293 subjects, 2711 sessions) to quantify microstructural features and length/volume features of several SWM systems. We find that axial, radial, and mean diffusivities show positive associations with age, while fractional anisotropy has negative associations with age in SWM throughout the entire brain. These associations were most pronounced in the frontal, temporal, and temporoparietal regions. Moreover, measures of SWM volume and length decrease with age in a heterogenous manner across the brain, with different rates of change in inter-gyri and intra-gyri SWM, and at slower rates than well-studied long-range white matter pathways. These features, and their variations with age, provide the background for characterizing normal aging, and, in combination with larger association pathways and gray matter microstructural features, may provide insight into fundamental mechanisms associated with aging and cognition.

Whole-brain, gray and white matter time-locked functional signal changes with simple tasks and model-free analysis.

Recent studies have revealed the production of time-locked blood oxygenation-level dependent (BOLD) functional MRI (fMRI) signals throughout the entire brain in response to a task, challenging the idea of sparse and localized brain functions, and highlighting the pervasiveness of potential false negative fMRI findings. In these studies, 'whole-brain' refers to gray matter regions only, which is the only tissue traditionally studied with fMRI. However, recent reports have also demonstrated reliable detection and analyses of BOLD signals in white matter which have been largely ignored in previous reports. Here, using model-free analysis and simple tasks, we investigate BOLD signal changes in both white and gray matters. We aimed to evaluate whether white matter also displays time-locked BOLD signals across all structural pathways in response to a stimulus. We find that both white and gray matter show time-locked activations across the whole-brain, with a majority of both tissue types showing statistically significant signal changes for all task stimuli investigated. We observed a wide range of signal responses to tasks, with different regions showing very different BOLD signal changes to the same task. Moreover, we find that each region may display different BOLD responses to different stimuli. Overall, we present compelling evidence that the whole brain, including both white and gray matter, show time-locked activation to multiple stimuli, not only challenging the idea of sparse functional localization, but also the prevailing wisdom of treating white matter BOLD signals as artefacts to be removed.

Anomalous and heterogeneous characteristics of the BOLD hemodynamic response function in white matter.

Detailed knowledge of the BOLD hemodynamic response function (HRF) is crucial for accurate analyses and interpretation of functional MRI data. Considerable efforts have been made to characterize the HRF in gray matter (GM), but much less attention has been paid to BOLD effects in white matter (WM). However, several recent reports have demonstrated reliable detection and analyses of WM BOLD signals both after stimulation and in a resting state. WM and GM differ in composition, energy requirements, and blood flow, so their neurovascular couplings also may well be different. We aimed to derive a comprehensive characterization of the HRF in WM across a population, including accurate measurements of its shape and its variation along and between WM pathways, using resting-state fMRI acquisitions. Our results show that the HRF is significantly different between WM and GM. Features of the HRF, such as a prominent initial dip, show strong relationships with features of the tissue microstructure derived from diffusion imaging, and these relationships differ between WM and GM, consistent with BOLD signal fluctuations reflecting different energy demands and neurovascular couplings in tissues of different composition and function. We also show that the HRF varies in shape significantly along WM pathways and is different between different WM pathways, suggesting the temporal evolution of BOLD signals after an event vary in different parts of the WM. These features of the HRF in WM are especially relevant for interpretation of the biophysical basis of BOLD effects in WM.

Joint independent component analysis for hypothesizing spatiotemporal relationships between longitudinal gray and white matter changes in preclinical Alzheimer's disease.

Characterizing relationships between gray matter (GM) and white matter (WM) in early Alzheimer's disease (AD) would improve understanding of how and when AD impacts the brain. However, modeling these relationships across brain regions and longitudinally remains a challenge. Thus, we propose extending joint independent component analysis (jICA) into spatiotemporal modeling of regional cortical thickness and WM bundle volumes leveraging multimodal MRI. We jointly characterize these GM and WM features in a normal aging (n=316) and an age- and sex-matched preclinical AD cohort (n=81) at each of two imaging sessions spaced three years apart, training on the normal aging population in cross-validation and interrogating the preclinical AD cohort. We find this joint model identifies reproducible, longitudinal changes in GM and WM between the two imaging sessions and that these changes are associated with preclinical AD and are plausible considering the literature. We compare this joint model to two focused models: (1) GM features at the first session and WM at the second and (2) vice versa. The joint model identifies components that correlate poorly with those from the focused models, suggesting the different models resolve different patterns. We find the strength of association with preclinical AD is improved in the GM to WM model, which supports the hypothesis that medial temporal and frontal thinning precedes volume loss in the uncinate fasciculus and inferior anterior-posterior association fibers. These results suggest that jICA effectively generates spatiotemporal hypotheses about GM and WM in preclinical AD, especially when specific intermodality relationships are considered a priori.

Multimodal neuroimaging in pediatric type 1 diabetes: a pilot multisite feasibility study of acquisition quality, motion, and variability.

Type 1 diabetes (T1D) affects over 200,000 children and is associated with an increased risk of cognitive dysfunction. Prior imaging studies suggest the neurological changes underlying this risk are multifactorial, including macrostructural, microstructural, and inflammatory changes. However, these studies have yet to be integrated, limiting investigation into how these phenomena interact. To better understand these complex mechanisms of brain injury, a well-powered, prospective, multisite, and multimodal neuroimaging study is needed. We take the first step in accomplishing this with a preliminary characterization of multisite, multimodal MRI quality, motion, and variability in pediatric T1D. We acquire structural T1 weighted (T1w) MRI, diffusion tensor MRI (DTI), functional MRI (fMRI), and magnetic resonance spectroscopy (MRS) of 5-7 participants from each of two sites. First, we assess the contrast-to-noise ratio of the T1w MRI and find no differences between sites. Second, we characterize intervolume motion in DTI and fMRI and find it to be on the subvoxel level. Third, we investigate variability in regional gray matter volumes and local gyrification indices, bundle-wise DTI microstructural measures, and N-acetylaspartate to creatine ratios. We find the T1-based measures to be comparable between sites before harmonization and the DTI and MRS-based measures to be comparable after. We find a 5-15% coefficient of variation for most measures, suggesting ~150-200 participants per group on average are needed to detect a 5% difference across these modalities at 0.9 power. We conclude that multisite, multimodal neuroimaging of pediatric T1D is feasible with low motion artifact after harmonization of DTI and MRS.

Minimal number of sampling directions for robust measures of the spherical mean diffusion weighted signal: Effects of sampling directions, b-value, signal-to-noise ratio, hardware, and fitting strategy.

Several recent multi-compartment diffusion MRI investigations and modeling strategies have utilized the orientationally-averaged, or spherical mean, diffusion-weighted signal to study tissue microstructure of the central nervous system. Most experimental designs sample a large number of diffusion weighted directions in order to calculate the spherical mean signal, however, sampling a subset of these directions may increase scanning efficiency and enable either a decrease in scan time or the ability to sample more diffusion weightings. Here, we aim to determine the minimum number of gradient directions needed for a robust measurement of the spherical mean signal. We used computer simulations to characterize the variation of the measured spherical mean signal as a function of the number of gradient directions, while also investigating the effects of diffusion weighting (b-value), signal-to-noise ratio (SNR), available hardware, and spherical mean fitting strategy. We then utilize empirically acquired data in the brain and spinal cord to validate simulations, showing experimental results are in good agreement with simulations. We summarize these results by providing an intuitive lookup table to facilitate the determination of the minimal number of sampling directions needed for robust spherical mean measurements, and give recommendations based on SNR and experimental conditions.

Bench to bore ramifications of inter-subject head differences on RF shimming and specific absorption rates at 7T.

PURPOSE: This study shows how inter-subject variation over a dataset of 72 head models results in specific absorption rate (SAR) and B1+ field homogeneity differences using common shim scenarios. METHODS: MR-CT datasets were used to segment 71 head models into 10 tissue compartments. These head models were affixed to the shoulders and neck of the virtual family Duke model and placed within an 8 channel transmit surface-loop array to simulate the electromagnetic fields of a 7T imaging experiment. Radio frequency (RF) shimming using the Gerchberg-Saxton algorithm and Circularly Polarized shim weights over the entire brain and select slices of each model was simulated. Various SAR metrics and B1+ maps were calculated to demonstrate the contribution of head variation to transmit inhomogeneity and SAR variability. RESULTS: With varying head geometries the loading for each transmit loop changes as evidenced by changes in S-parameters. The varying shim conditions and head geometries are shown to affect excitation uniformity, spatial distributions of local SAR, and SAR averaging over different pulse sequences. The Gerchberg-Saxton RF shimming algorithm outperforms circularly polarized shimming for all head models. Peak local SAR within the coil most often occurs nearest the coil on the periphery of the body. Shim conditions vary the spatial distribution of SAR. CONCLUSION: The work gives further support to the need for fast and more subject specific SAR calculations to maintain safety. Local SAR10g is shown to vary spatially given shim conditions, subject geometry and composition, and position within the coil.

Detection of functional activity in brain white matter using fiber architecture informed synchrony mapping.

A general linear model is widely used for analyzing fMRI data, in which the blood oxygenation-level dependent (BOLD) signals in gray matter (GM) evoked in response to neural stimulation are modeled by convolving the time course of the expected neural activity with a canonical hemodynamic response function (HRF) obtained a priori. The maps of brain activity produced reflect the magnitude of local BOLD responses. However, detecting BOLD signals in white matter (WM) is more challenging as the BOLD signals are weaker and the HRF is different, and may vary more across the brain. Here we propose a model-free approach to detect changes in BOLD signals in WM by measuring task-evoked increases of BOLD signal synchrony in WM fibers. The proposed approach relies on a simple assumption that, in response to a functional task, BOLD signals in relevant fibers are modulated by stimulus-evoked neural activity and thereby show greater synchrony than when measured in a resting state, even if their magnitudes do not change substantially. This approach is implemented in two technical stages. First, for each voxel a fiber-architecture-informed spatial window is created with orientation distribution functions constructed from diffusion imaging data. This provides the basis for defining neighborhoods in WM that share similar local fiber architectures. Second, a modified principal component analysis (PCA) is used to estimate the synchrony of BOLD signals in each spatial window. The proposed approach is validated using a 3T fMRI dataset from the Human Connectome Project (HCP) at a group level. The results demonstrate that neural activity can be reliably detected as increases in fMRI signal synchrony within WM fibers that are engaged in a task with high sensitivities and reproducibility.

Aging and white matter microstructure and macrostructure: a longitudinal multi-site diffusion MRI study of 1218 participants.

Quantifying the microstructural and macrostructural geometrical features of the human brain's connections is necessary for understanding normal aging and disease. Here, we examine brain white matter diffusion magnetic resonance imaging data from one cross-sectional and two longitudinal data sets totaling in 1218 subjects and 2459 sessions of people aged 50-97 years. Data was drawn from well-established cohorts, including the Baltimore Longitudinal Study of Aging data set, Cambridge Centre for Ageing Neuroscience data set, and the Vanderbilt Memory & Aging Project. Quantifying 4 microstructural features and, for the first time, 11 macrostructure-based features of volume, area, and length across 120 white matter pathways, we apply linear mixed effect modeling to investigate changes in pathway-specific features over time, and document large age associations within white matter. Conventional diffusion tensor microstructure indices are the most age-sensitive measures, with positive age associations for diffusivities and negative age associations with anisotropies, with similar patterns observed across all pathways. Similarly, pathway shape measures also change with age, with negative age associations for most length, surface area, and volume-based features. A particularly novel finding of this study is that while trends were homogeneous throughout the brain for microstructure features, macrostructural features demonstrated heterogeneity across pathways, whereby several projection, thalamic, and commissural tracts exhibited more decline with age compared to association and limbic tracts. The findings from this large-scale study provide a comprehensive overview of the age-related decline in white matter and demonstrate that macrostructural features may be more sensitive to heterogeneous white matter decline. Therefore, leveraging macrostructural features may be useful for studying aging and could facilitate comparisons in a variety of diseases or abnormal conditions.

Functional Parcellation of Human Brain Using Localized Topo-Connectivity Mapping.

The analysis of connectivity between parcellated regions of cortex provides insights into the functional architecture of the brain at a systems level. However, the derivation of functional structures from voxel-wise analyses at finer scales remains a challenge. We propose a novel method, called localized topo-connectivity mapping with singular-value-decomposition-informed filtering (or filtered LTM), to identify and characterize voxel-wise functional structures in the human brain from resting-state fMRI data. Here we describe its mathematical formulation and provide a proof-of-concept using simulated data that allow an intuitive interpretation of the results of filtered LTM. The algorithm has also been applied to 7T fMRI data acquired as part of the Human Connectome Project to generate group-average LTM images. Generally, most of the functional structures revealed by LTM images agree in the boundaries with anatomical structures identified by T1-weighted images and fractional anisotropy maps derived from diffusion MRI. In addition, the LTM images also reveal subtle functional variations that are not apparent in the anatomical structures. To assess the performance of LTM images, the subcortical region and occipital white matter were separately parcellated. Statistical tests were performed to demonstrate that the synchronies of fMRI signals in LTM-derived functional parcels are significantly larger than those with geometric perturbations. Overall, the filtered LTM approach can serve as a tool to investigate the functional organization of the brain at the scale of individual voxels as measured in fMRI.

Dynamic variations of resting-state BOLD signal spectra in white matter.

Recent studies have demonstrated that the mathematical model used for analyzing and interpreting fMRI data in gray matter (GM) is inappropriate for detecting or describing blood-oxygenation-level-dependent (BOLD) signals in white matter (WM). In particular the hemodynamic response function (HRF) which serves as the regressor in general linear models is different in WM compared to GM. We recently reported measurements of the frequency contents of resting-state signal time courses in WM that showed distinct power spectra which depended on local structural-vascular-functional associations. In addition, multiple studies of GM have revealed how functional connectivity between regions, as measured by the correlation between BOLD time series, varies dynamically over time. We therefore investigated whether and how BOLD signals from WM in a resting state varied over time. We measured voxel-wise spectrograms, which reflect the time-varying spectral patterns of WM time courses. The results suggest that the spectral patterns are non-stationary but could be categorized into five modes that recurred over time. These modes showed distinct spatial distributions of their occurrences and durations, and the distributions were highly consistent across individuals. In addition, one of the modes exhibited a strong coupling of its occurrence between GM and WM across individuals, and two communities of WM voxels were identified according to the hierarchical structures of transitions among modes. Moreover, these modes are coupled to the shape of instantaneous HRFs. Our findings extend previous studies and reveal the non-stationary nature of spectral patterns of BOLD signals over time, providing a spatial-temporal-frequency characterization of resting-state signals in WM.